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The refractive index and electronic gap of water 
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Determining the electronic and dielectric properties of water at high pressure and 
temperature is an essential prerequisite to understand the physical and chemical properties 
of aqueous environments under supercritical conditions, for example, in the Earth interior. 
However, optical measurements of compressed ice and water remain challenging, and it has 
been common practice to assume that their band gap is inversely correlated with the 
measured refractive index, consistent with observations reported for hundreds of materials. 
Here we report ob initio molecular dynamics and electronic structure calculations showing 
that both the refractive index and the electronic gap of water and ice increase with increasing 
pressure, at least up to 30GPa. Subtle electronic effects, related to the nature of interband 
transitions and band edge localization under pressure, are responsible for this apparently 
anomalous behaviour. 
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Water is arguably the most important material in the 
biosphere, as well as a key constituent of the Earth 
crust and mantle. Its properties have been extensively 
studied as a function of temperature (T) and pressure (P) for 
several decades 1 ' 2 . However, the electronic structure of water 
under pressure has remained elusive, due to experimental 
difficulties in measuring absorption processes in wide band gap 
insulators in diamond anvil cells. The electronic properties of 
compressed water are key, for example, to understanding electron 
charge transfer rates in redox reactions 3 ' 4 occurring in 
supercritical water 1 and aqueous solutions in the Earth's 
interior 2 ; such reactions, which depend on both the electronic 
and static dielectric constant of water, determine the oxidation 
states of minerals and rocks. 

Water was predicted to undergo an insulator-to-metal transi- 
tion at rather large pressures, estimated between 100 GPa and 
5TPa 5-8 , although no experimental confirmation has yet been 
reported. Experimentally, it is difficult to directly establish the 
variation of the liquid electronic gap as a function of P, which at 
ambient conditions is outside the optical window of diamond: the 
optical gap of diamond is 5.4 eV (ref. 9) with the absorption tail as 
low as ~4eV, and the quasi particle gap of water at ambient 
conditions is estimated to be 8.7±0.5eV (refs 10,11) with the 
optical absorption onset at ~7eV (refs 12,13). Hence, the band 
gap (Eg) of high pressure ice (a disordered ice phase stable 
between 3 and ~ 60 GPa (refs 14,15)) has been inferred from 
measurements of its refractive index, n, which can be obtained by 
using low energy photons without probing absorption 
processes 16 . In particular, it has been common practice to use a 
single- oscillator model, for example, the Penn model 17 , to relate n 
and Eg, as in the case of other molecular crystals, notably 
hydrogen 18 ' 19 . According to the Penn model, the refractive index 
n and the electronic gap E g are inversely correlated: 

n 2 = 1 + where co fp is the plasma frequency. This 

inverse correlation between n and E g holds for a broad class of 
materials 20 ' 21 . 

Zha et al. 16 measured the refractive index of ice VII as a 
function of pressure up to 120 GPa and found that it increased 
with pressure, similar to the case of solid hydrogen. They found 
weak dispersion of the refractive index, and its pressure 
dependence, leading them to conclude that the band gap is 
preserved to very high pressure, ultimately closing only in the 
terapascal pressures and well beyond the stability field of ice VII. 
The refractive index of liquid water has been measured only up to 
7 GPa, at T = 673 K (ref. 22), while no measurements have yet 
been reported of its electronic gap in such P-T regime. 

Here we report ab initio molecular dynamics (MD) simulations 
of water under pressure and calculations of the refractive index n 
and the electronic gap E g in the pressure range 0-30 GPa. The 
values obtained for the refractive index of ice and water 
are in good agreement with experiments 16 ' 22 ; however, we find 
that both n and E g increase with increasing pressure in the solid 
and liquid phases. Our results show that a single-oscillator model 
used to rationalize the electronic properties of many 
semiconductors and insulators may not be employed to 
describe water under pressure. We find that, because of subtle 
but important changes in the localization properties of the 
valence band edge of water and ice, the strength of inter band 
transitions depends on pressure. As a consequence n and E g are 
not inversely correlated, but exhibit a more complex 
interdependence, not captured by widely used simple models. 
Our results are consistent with those of Hermann and 
Schwerdtfeger 23 who reported an increase of the optical gap of 
ice under pressure, using calculations based on many body 
perturbation theory with a fitted dielectric constant. 
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Results 

Refractive index of water and ice. For photon energies (no) much 
smaller than the band gap, absorption processes may be ignored and 
n — y/ei, where c Y is the real part of the electronic dielectric con- 
stant. For simple semiconductors, the Penn model 17 yields a 
reasonable approximation of ex. The validity of such model may 
extend beyond simple semiconductors, as suggested by Wemple and 
DiDomenico 20 who showed that [ei(0) — 1] is approximately 
proportional to l/E g for more than 100 solids and liquids. 
However, as discussed below, the description of the properties of 
water and ice requires a more sophisticated level of theory. 

Figure 1 reports the calculated refractive index of water and ice, 
as a function of pressure, and shows that n increases with 
increasing P along isotherms for both the solid and the liquid. We 
also show experimental results for ice VII (ref. 16) and water up 
to 7 GPa (ref. 22). To the best of our knowledge, our results 
represent the first data obtained for the refractive index of the 
supercritical liquid above 7 GPa. Indeed, the highly corrosive 
character of hot, compressed water makes experiments rather 
challenging at high P and T. 

We find good agreement with experiments, especially above 
5 GPa, with the deviations present at low pressure ascribed to 
errors introduced by the use of the Perdew-Burke-Ernzerhof 
(PBE) functional 24 . Semi-local functionals have been shown to 
yield results for water and ice under pressure in better agreement 
with experiments than at ambient conditions 25 , for example, for 
the calculation of the equation of state 26 and the dielectric 
properties 15 ' 26 . Preliminary calculations of the electronic 
dielectric constant of water at ambient conditions using the 
PBE0 (refs 27,28) functional showed a much improved agreement 
with experiment, yielding e 1 = 1.78 (to be compared with the 
experimental value of 1.77). The same improvement was also 
found for the static dielectric constant of ice Ih 29 . For ice VIII at 
30 GPa, we find instead a minor difference of 4% between the 
electronic dielectric constants obtained by PBE and PBE0. 
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Figure 1 | Refractive indices. Pressure dependence of the refractive index, 
n, of water and ice. We computed n = where e y is the real part of the 
electronic dielectric constant, obtained from ab initio simulations (see text). 
The experimental data for ice VII, obtained with light wavelength of 
630 nm, are from ref. 16. The experimental values for water at 673 K are 
from ref. 22. We note that the statistical fluctuations in the value of n are 
rather small, indicating that e\ is largely insensitive to reorientations 
of the molecules in the fluid and to hydrogen bonding rearrangements, 
unlike the ionic dielectric constant 26 . When increasing the temperature 
from 1,000 to 2,000 K, n showed again moderate variations. The error bars 
show the s.d. values of the refractive indices of water. Lines are 
drawn to guide the eyes only. Calculations were conducted for ice VIII 
instead of VII for computational simplicity, as ice VII is proton disordered. 
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Electronic gap of water and ice. In Fig. 2, we show the electronic 
gap of water and ice, computed as the difference between the single 
particle energies of the conduction band minimum (CBM) and the 
valence band maximum (VBM) as a function of pressure. We 
report values obtained with the PBE functional and, in the inset, 
with the hybrid functional, PBEO (refs 27,28). In spite of large 
statistical fluctuations at high T, the trend of an increasing gap with 
pressure is evident, irrespective of the level of theory used. As well 
known, the PBE functional underestimates the band gap of 
water 11 ' 30 , but the gap variation under pressure is qualitatively 
the same and quantitatively similar within PBE and PBEO. We 
expect the rate at which the optical gap increases with pressure be 
larger than for the electronic gap. Indeed, based on our results for 
£1, we estimated that the exciton binding energy, £ b , of water will be 
decreased by about 60% in going from ambient pressure to 
-30 GPa, and that of ice VII/VIII by 40% (The E h of water at 
ambient conditions is about 2.4 eV (ref. 30)). Our estimate is based 
on the relation E^ — m* / (m e c\ao) , where ra* is the electronic 
effective mass, ra e the electron mass and a 0 the Bohr radius 31 , and 
we assumed a negligible variation of the electronic effective mass 
under pressure, as indicated by our band structure calculations (not 
shown). Our results for band gaps are consistent with those 
reported by Hermann and Schwerdtfeger 23 for ice VIII, using many 
body perturbation theory. These authors used the G 0 W 0 
approximation and the Bethe-Salpeter equation 23 to compute 
absorption spectra, with a model dielectric function to approximate 
the screened Coulomb interaction. In their model, the electronic 
dielectric constant was an input parameter. Our results are also 
consistent with those of Boero et al 32 who reported an increase in 
the band gap of supercritical water at 653 K, for densities larger than 
0.5 g cm _ 3 , by using ab initio MD and the BLYP functional. 

Relation between refractive index and band gap. Our calcula- 
tions showed that both the refractive index and the electronic gap of 
water and ice increase under pressure, up to at least 30 GPa. They 
also showed that n of ice is larger than that of the liquid, in spite of 
the latter having a smaller electronic gap. Thus, the simple inverse 
correlation between E g and n used to interpret experiments, and valid 
for a wide range of substances 20 , does not hold for water and ice. We 
show in the following that two reasons are responsible for this 
apparent anomalous behaviour: an increase in the electronic density 



with pressure and hence of the plasma frequency cd\ p = An ^ e , where 
p e is the density of valence electrons and e the elementary charge; a 
change in the localization property of the valence band of the liquid 
and solid, as P is increased, which in turn is responsible for changes 
in the strength of interband transitions. 

Interestingly, while for ice the increase in p e (and C0fp) 
counterbalances the increase in E g and hence [£i(0) — l]/p e 
decreases under pressure, the corresponding quantity for water 
shows an increase (see Fig. 3). Ice and water behave as several 
oxides 33 in exhibiting a positive derivative of the refractive index, 
with respect to p e . However, for both water and ice the quantity 
[ei(0) — l]Eg/pg is not a constant, as it would be if the Penn 
model correctly described the relation between the dielectric 
constant and electronic gap, but it monotonically increases with 
pressure, at all T considered here. 

We therefore adopted a higher level of theory and derived the 
real part of the electronic dielectric constant in the random phase 
approximation: 

61 (CO) = 1 + -4-Y [ dk ^ (1) 
' BZ CV 

where c and v are the conduction and valence band indices, 
respectively, and the integral is over the Brillouin Zone (BZ). In 
the integrand,/^ (k J is the oscillator strength in the polarization 
direction associated with the transition between the bands c 
and v, and h co CY (k) is the corresponding transition energy. In 
the case of water, only the T point was used in our electronic 
structure calculations; at zero frequency one has: 

£l (o) = i + ^y4 (2) 

where V is the volume of the simulation cell. Note that in a single- 
oscillator model with tico cv = £ g , Equation 2 reduces to the Penn 
model, once the oscillator strength sum rule is applied. 
The oscillator strength is given by 

^ = ^^ im . 7]D)] 2 x£2{(0er)Vo)er> (3) 

where <c| and |i)> denote Kohn-Sham orbitals corresponding to 
valence (v) and conduction (c) states, and e 2 is the imaginary part 
of the dielectric function. For v to c transitions right across the 
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Figure 2 | Electronic gaps. Computed electronic gap of water and ice VIII 
as a function of pressure at 0, 1,000 and 2,000 K. The calculations were 
performed using the semilocal density functional, PBE, and in the inset, the 
hybrid functional PBEO. The error bars show the s.d. values of the band gaps 
of water. Lines are drawn to guide the eyes only. 
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Figure 3 | Volume effect on refractive indices. This figure shows G = 
(ei(0)-1)/p e as a function of pressure at 0, 1,000 and 2,000 K, where 
e-i(O) is the real part of the electronic dielectric constant at zero frequency 
and p e the valence electron density. The error bars show the s.d. values of 
the data for water. Lines are drawn to guide the eyes only. 
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band gap, the oscillator strength is proportional to [^(gw) Vcd^]. 
In Fig. 4 we plotted e 2 (co)V of two water configurations at 
1,000 K, at ~ 1 GPa and ~ lOGPa. These are representative of a 
set of 25 and 55 snapshots analysed for each pressure, extracted 
from our MD trajectories at equally spaced simulation times; for 
each of them we found that at nco = £ g , e 2 V is much larger at 
~ 10 GPa than at ~ 1 GPa. The oscillator strength of the 
transitions just across the band gap, which yields the largest 
contribution to e^O) according to Equation 2, is thus substantially 
enhanced when the pressure increases from ~ 1 GPa to ~ 10 GPa, 
and the absorption edge is blue shifted, as expected from our 
electronic gap calculations. We conclude that it is because of the 
enhancement off£, which outweighs the increase of £ g , that the 
refractive index of water increases under pressure. 

We note that local field effects are not included in Equation 1 
(ref. 34); they are instead taken into account in our density 
functional perturbation theory (DFPT) 35 calculation of the 
refractive index (Denote the matrix elements of the static 
dielectric matrix e GG >(0), where G and G' are reciprocal lattice 
vectors; ef 1 (0) = [e GG < (0)] ~ 1 1 G=G , =0 > implying that ef 1 (0) is the 
first element of the inverse dielectric matrix (ref. 35)). However, 
in Equation 1, e^O) is simply e 00 (0) and so called local field effects 
are neglected 34 ). We verified that neglecting such effects amounts 
to a negligible error in the calculations of n. In Fig. 5, we compare 
£i(0) of ice VIII, as well as that of high pressure water, obtained 
with and without local field corrections. The two sets of values are 
rigidly shifted by only 0.1 ~ 0.2 with respect to each other, 
indicating that none of the trends as a function of pressure 
reported here is affected by the inclusion of local field effects. 

Localization properties of water and ice valence band edges. To 

understand the changes in the oscillator strength of water under 
high pressure, we computed the inverse participation ratios 
(IPRs) of the VBM and CBM orbitals; they are reported in Fig. 6. 
The IPR is defined as 



IPRa 



(4) 



where N is the total number of points used to perform integrals 
over a real space grid, and i// a (i) is the wave function associated 




Figure 4 | Oscillator strength of interband transitions. The imaginary part 
of the electronic dielectric constant of water, e 2 multiplied by the cell 
volume Vas a function of frequency co at ~1 and ~10 GPa and 1,000 K. E g 
labels the frequencies corresponding to the band gaps of water at ~1 and 
~10GPa, respectively. 



with the band oc. The quantity IPR a yields the statistical variance 
for the distribution of \i// a | 2 . The larger the IPR, the more loca- 
lized the wave function. An IPR a close to 1 indicates a delocalized 
electronic orbital. Figure 6 shows that, under pressure, as at 
ambient conditions, the CBM of water is much more delocalized 
than its VBM 37 . The interesting finding is that with increasing 
pressure the localization properties of the CBM are largely 
unaffected, while those of the VBM change substantially, with a 
tendency of the VBM to delocalize. The electronic states of ice 
VIII behave similarly except that the variation of the VBM 
localization as a function of P is smaller than that of water, as 
shown in Fig. 6c. As a result, under pressure the spatial overlap 
between the electron and hole wave functions is enhanced, and 
the matrix element (c\e£ • 7*1 d) increases (see Eq. 3). The fact 
that the localization of the CBM does not change significantly 
under pressure, while that of the VBM substantially changes, 
suggests that the opening of the gap in ice and water under 
pressure might be mostly due to a change in the position of the 
VBM. 

Experimentally Zha et al. 16 measured n as a function of co and 
obtained two parameters, E d and E 0 from the equation n 2 — 1 = 
E d E 0 /(Eo — h 2 co 2 ). The fit was then used to conclude that the gap 
of ice decreases under pressure. Note that in the experiment of 
ref. 16 when the light wavelength is between 569 to 741 nm, n 
changes by < 0.05; as a result, the effective band gap determined 
from the single-oscillator model changes little with pressure. It is 
presumably difficult to estimate E d and E 0 . The values for E d and 
E 0 obtained from the fits gave a predicted closing of the band gap 
only at ultrahigh pressures ( ~4.5 TPa), well beyond the pressure 
range of the measurements. We have shown here that state-of- 
the-art calculations of the electronic structure are essential to 
properly interpret complex and sophisticated optical experiments 
and obtain the correct band gap trend under pressure. 



Discussion 

We reported the first data for the refractive index of supercritical 
water above 7 GPa and at high temperature. We showed that the 
refractive index, and the electronic and optical gaps of water and ice 
increase under pressure, and thus water properties differ 
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Figure 5 | Local field effects on dielectric constants. The dashed lines 
show the real part of the electronic dielectric constants of water and ice VIII 
(0 K) at zero frequency, ^(0), obtained from density functional 
perturbation theory (DFPT) without local field effects (NLF). As a 
comparison, the results obtained with DFPTand including local field effects 
(as done in Fig. 1), are shown by the solid lines. The error bars show the s.d. 
values of the electronic dielectric constants of water. 
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Figure 6 | Localization properties of water conduction and valence band edges. The inverse participation ratios (I PR) of the valence band maximum 
(VBM, circles) and conduction band minimum (CBM, squares) of (a,b) water and (c,d) ice VIII as a function of pressure. Results of panels b and d were 
obtained with the hybrid functional PBE0, while those of panels a and c with the semi-local functional PBE. The error bars show the s.d. values of the IPRs of 
water. Lines are drawn to guide the eyes only. 



substantially from those of apolar molecular fluids and solids, for 
example, hydrogen 38 , methane 39,40 and benzene 41,42 , whose 
electronic gaps decrease with pressure but refractive indices 
increase. 

We showed that the unusual relationship between n and E g 
stems from an enhancement of the strength of interband 
transitions right across the band gap. This enhancement, in 
combination with a decrease of the volume, outweighs the 
increase in the band gap upon compression, and hence the 
refractive index of water increases with increasing pressure. The 
behaviour of the interband transition strength originates from the 
derealization of the valence band edge, whose overlap with the 
CBM increases with pressure. Our results showed that the use of 
simple models, for example, the Penn model 17 , to relate measured 
refractive indices to band gaps of molecular fluids and crystals 
should be revisited and in general complex electronic structure 
effects should be taken into account, when predicting values of 
band gaps. We emphasize that the type of electronic structure 
calculations reported here allows one to infer the band gap of 
water from high pressure measurements of the refractive index 16 . 
Measurements of the band gap of water by light absorption or 
emission are not yet feasible under pressure, in diamond anvil 
cells; therefore, establishing the correlation between the band gap 
and the refractive index is a crucial step to infer band gaps from 
measured dielectric constants. In general, ab initio calculations 
play a key role in interpreting measurements of fluids under 
pressure, especially in unraveling their electronic structure. 

Finally, we note that the results reported here for gaps and 
dielectric constants are expected to have an important impact in 
the broad field of water-related science. For example, factors 
determining the rate of electron transfer in aqueous environments 
(and hence several chemical reactions occurring in these media) 
include the refractive index and the static dielectric constant of 



the medium 3 . Hence, knowledge of such quantities is key to 
predict reorganization energies and electron transfer rates in 
supercritical water and water under pressure, for example, in the 
Earth's interior. In addition, the mobility of charges created in 
water by, e.g., photo -excitation, depends, among other factors, on 
the exciton binding energy (£ b ); our results indicate that E h 
decreases under pressure and hence charge mobility is expected to 
change substantially in supercritical water, with respect to water 
at ambient conditions, with implications for charge transport 
under bias, in supercritical conditions. 

Methods 

Water. We carried out ab initio MD simulations of water under pressure as a 
function of temperature in the Born-Oppenheimer approximation with the Qbox 
code (version 1.54.2, http://eslab.ucdavis.edu/software/qbox/) 43 . We used 128- 
molecule unit cells with periodic boundary conditions, plane-wave basis sets and 
norm conserving pseudopotentials (Gygi, F. Pseudopotential Table, http:// 
fpmd.ucdavis.edu/potentials/index.htm). Atomic trajectories were generated with a 
plane-wave kinetic energy cutoff of 85 Ry, which was then increased to 220 Ry for 
the calculation of pressure on selected snapshots 26 . 

Ice VIII. Calculations for ice VIII were conducted by using a 8-molecule unit cell 
with a plane-wave energy cutoff of 220 Ry and 4x4x3 k-point mesh and a 96- 
molecule supercell with a cutoff of 85 Ry and T point only. 

Electronic structure calculations. Trajectories were generated with the PBE 
exchange correlation functional 24 , while electronic gaps and dielectric constants 
were computed with both the PBE and PBE0 (refs 27,28) functionals on snapshots 
extracted from PBE trajectories. At the PBE level, the electronic dielectric constant 
Ei was calculated using density functional perturbation theory (DFPT) 36 , as 
recently implemented in the Qbox code for the calculation of Raman spectra 44 . 
At the PBE0 level of theory, we instead used an electronic enthalpy functional 45 , 
where the polarization was computed using maximally localized Wannier 
functions 46 with a correction scheme 47 . We calculated the strengths of inter band 
transitions with the Quantum-expresso package 48 . We used a Gaussian broadening 
parameter equal to 0.001 eV. 
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